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Wilson's momentum shell renormalization group method is used to solve for the dynamics of the 
dissipative two-state system. We utilize the mapping of the spin-boson model onto the anisotropic 
Kondo model (AKM) and solve for the dynamics of the latter. We find that the AKM captures 
the physics of the dissipative two-state system for dissipation strength < a < 4 corresponding to 
°° ^ -^ll ^ — °° ™ the AKM. The dynamics of the AKM shows a smooth crossover between two 
strong coupling regimes corresponding to weak and strong dissipation in the spin-boson model. 
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A problem which is of general interest in both physics 
and chemistry is that of a quantum mechanical system 
tunneling between two states and subject to a dissipative 
coupling to environmental degrees of freedom. Examples 
of this type of system abound; they include the tunnel- 
ing of defects in metallic glasses [1], the motion of the 
total flux in a SQUID between two metastable fluxoid 
states [2] and the diffusion of protons and muons in met- 
als [3]. The main theoretical interest lies in a description 
of the dynamics of the generalized coordinate of the two 
level system subject to the influence of the environment. 
Both the equilibrium and non-equilibrium dynamics are 
of interest for the different experimental realizations of 
two-level systems. In the case of macroscopic quantum 
coherence experiments in SQUIDs, the system can be 
prepared in one of the two states by applying a bias ( an 
external magnetic field) for times i < and then allowed 
to evolve for i > in zero bias. The non-equilibrium 
correlation functions of the two state system are then of 
primary interest. For most microscopic systems such an 
initial state preparation is not realizable and the interest 
then lies in the equilibrium dynamics. In this paper we 
present results for the equilibrium case with a bias. 

Specifically we consider the spin-boson Hamiltonian, 

HsB = -^hAa^ + ^ea^ + ^^^(a^ac + ^) 

a 

+ 0^0^' J2 /O " («a +4)- (1) 

a 

Here cr^, i = x,y, z are Pauli spin matrices, the two-states 
of the system correspond to =| and =|. A is the 
bare tunneling matrix element and e is a bias. The en- 
vironment is represented by an infinite set of harmonic 
oscillators (labelled by the index a) with masses rria and 
frequency spectrum coupling linearly to the coordi- 
nate Q = ^qoCTz of the two-level system via a term char- 
acterized by the couplings C^. The environment spectral 
function is given in terms of these couplings, oscillator 



iVa). In the case of an Ohmic heat bath, of interest to us 
here, we have = 27raa;, for uj << w^, where is a 

high energy cut-off and a is a dimensionless parameter 
characterizing the strength of the dissipation. 

A great deal of work has been carried out over the 
last 10 years in order to understand the dynamics of this 
apparently simple model [4,5]. Extensive calculations 
based on the Feynman- Vernon path integral formal- 
ism within the "Non-Interacting Blip Approximation" 
(NIBA) have yielded reliable information for weak dissi- 
pation and short times [4] and have also provided some 
insight into the expected behaviour in other regimes. In 
addition to such direct attempts at calculating dynam- 
ics for the spin-boson model, it has proved fruitful to 
exploit analogies between this model and several other 
models, including the inverse-square Ising model [6], 
the anisotropic Kondo model (AKM) and the Vigmann- 
Finkel'stein (VF) model [7]. In the long time approxi- 
mation it has been shown that the partition functions, 
and hence the thermodynamics, of these different models 
can be put into correspondence and the parameters of 
the models related. This correspondence has also been 
carried out for the respective fermionic Hamiltonians by 
applying an approximate bosonization method valid for 
low energies uj « uic [8]. In this paper we exploit such 
a mapping of the spin-boson model onto the AKM in or- 
der to make predictions about the dynamics of the former 
on the basis of renormalization group calculations on the 
latter. It has not been clear [4] to what extent the above 
mapping is valid in all parameter regimes, in particular 
for weak dissipation, and this has partly motivated this 
work. The AKM is given by 
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masses and frequencies by J(w) = ■|X^(j( 



kk' 

where the first term represents non-interacting conduc- 
tion electrons and the second and third terms represent 
an exchange interaction between a localized spin 1/2 
and the conduction electrons with strength Jx, Jy. For 
Jx = J|| the model reduces to the usual Kondo model of 
magnetic impurities in metals. A local magnetic field, h, 
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coupling only to the impurity spin in the Kondo model 
(the last term in Eq. 2) corresponds to a finite bias, e, 
in the spin-boson model. The correspondence between 
H and Hsb is then given by e = gfish, ^ = pJ± and 

a = (1 + ^)'^, where tanS = — — ^-"-j 6 is the phase 
shift for scattering of electrons from a potential J||/4 and 
p = 1/2D is the conduction electron density of states per 
spin at the Fermi level for a flat band of width 2D [4]. We 
note that whereas the spin-boson model can describe dis- 
sipation with arbitrary strength, < a < co, the AKM 
is only capable of describing the region < a < 4 with 
a = corresponding to Jy = +co, a = 1 corresponding 
to the antiferromagnetic/ ferromagnetic boundary of the 
AKM (J|| = 0), and 1 < a < 4 corresponding to the 
ferromagnetic regime of the AKM. 

The relevant dynamical quantity for the two-level sys- 
tem is the response function xsb{^,T) = {{cr z] <^ z)) 
which translates into the local dynamic spin suscep- 
tibility, x('^i^) = {{Sz]Sz)), for the Kondo model. 
In order to calculate this quantity we apply Wilson's 
momentum shell renormalization group method which 
has recently been generalized to the calculation of dy- 
namical quantities for a number of models (e.g. [11]). 
Briefly, the procedure, explained in detail in [9], con- 
sists of (i) linearizing the spectrum about the Fermi 
energy ek — > vpk, (ii) introducing a logarithmic mesh 
of k points kn = A~" and (iii) performing a unitary 
transformation of the Cka such that foa = X^j. Cka is 
the first operator in a new basis, fna^n = 0,1,..., 
which tridiagonalizes He = Yjk^^k^c\^Ck^ in k-space, 

i.e. ^ Y.,Y.n=oU^-'''\fi+i,U + h.c), with 

^„ ^ (1 + A-i)/2 for n » 1. The Hamiltonian (2) 
with the above discretized form of the kinetic energy is 
now diagonalized by the following iterative process: (a) 
one defines a sequence of finite size Hamiltonians = 

E,El"o ^nA-"/2(/t^,^/„, + h.C.) + ^(/o\/oi5- + 

/oVoT^+) + ^(/ot/oT - foifoi)S^ for N > 0; (b) the 
Hamiltonians Hjv are rescaled by A = such that the 
energy spacing remains the same, i.e. Hjv = A = -ffjv- 
This defines a renormalization group transformation 
Hn+1 = A^/^Hjv + J2^L^N{flr+l^fN^L + h.c.) - Eg,n+i, 
with Eg^n+i chosen so that the ground state energy of 
Hn+1 is zero. Using this recurrence relation, the se- 
quence of Hamiltonians H^ for 7\r = 0, 1, . . . is iteratively 
diagonalized within a product basis of, typically, up to 
1200 states for each iteration. This gives the excitations 
and many body eigenstates at a corresponding set of en- 
ergy scales wjv defined by wjv = A = and allows a 
direct calculation of response functions. Thus, for exam- 
pie, x('^, T) is given by x{'^,T) = x'('^, T) + ix"{i^, T) = 

^Em,n \M^,n\^ ^+7o-{el-Z) ' "^"^"^ ^""^ many- 

body excitations of -ffjv, Zn{T) the corresponding par- 
tition function of -ffjv, and -M^„ =< mjS'^ln >jv" the 
relevant many-body matrix elements for the dynamic 



susceptibility. The quantity we actually calculate is 
S(uj) = —I X (^+"?) -yyjjjgjj jg related to the neutron scat- 
tering cross-section. Our results were obtained for A = 2, 
keeping the 320 lowest states at each iteration. In this 
paper we discuss only the T = results. 

The accuracy of the numerical calculations could be 
tested by (a) showing that the exactly solvable Toulouse 
limit, a = |, could be reproduced (described below) and 
(b) by verifying the generalized Shiba relation for the dy- 
namic spin susceptibility of the spin-boson model [12]. 
The latter relation states that, at T = 0, S[lli = 0) = 
2ax'(w = 0) . Generalized with the extra factor of a, it 
is also valid for the AKM, as explicitly verified by our 
numerical results. We chose to extract the static sus- 
ceptibility via a Kramer's-Kronig relation x'{^ = 0) = 
S[uj)duj. Verified in this way, the generalized Shiba 
relation provides a good test of the method, not just at 
low frequency, but at all energy scales and for arbitrary 
values of a and A. The error in x'C^) typically 5% 
and contributed the main source of error, approximately 
10%, to the Shiba relation; representative cases are shown 
in Table I [13]. In all cases we found that x"{^) ~ ^ for 
a; — > showing that the spin-spin correlation function, 
{[az{t), az{0)]), decays as for long times i.e. the tun- 
neling is always incoherent at long times. In contrast, the 
NIBA gives an algebraic decay ~ l/t'^^^~"') depending on 
the coupling a. 
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FIG. 1. The response function 5(cj) of the AKM for 
a = 0.1 (J|| = 4.698) and A = 0.08 (solid line), A = 0.1 
(dotted line), A = 0.15 (dashed line). The position of the 
inelastic peaks is at AJ! = 0.045, 0.063 and 0.087 respectively. 
The inset shows that 5(cj)/5(0) is a universal function of the 
rescaled energy w/A* (5(0) ~ -^irr)- Energies are in units of 
D=^ = l. 

The case of zero dissipation, a = 0, corresponds to 
J|| = CO in the AKM. This is exactly solvable; S[lli) = 
-^[6{uj — J±) — 6{uj + J±)] is a sum of two delta functions 
at a; = ±J± . In the spin-boson model the response is 
also a sum of two delta functions at a; = ±A. We can 
thus identify the cut-off appearing in — = pJ± as the 
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bandwidth, uic = 2D. In this case coherent oscillations 
with frequency A are realized. On increasing a, Fig. (1), 
the above peaks survive but are broadened. The sys- 
tem exhibits oscillations for short times t ^ 1/ , where 
A^ is the renormalized tunneling frequency (see below). 
As discussed above, the long time behaviour is always 
incoherent for any finite a and A. The position of the 
peaks, which is a measure of the renormalized tunnel- 
ing frequency, is reduced i.e. as expected, dissipation 
hinders tunneling. Fig. (1) illustrates these features for 
a = 0.1. For a << 1 and small A, scaling arguments 
give a renormalized tunneling frequency A^ = A(^) 
[4]. Numerical results for a = 0.01 gave for the inelastic 
peak positions. A*, values within 5% of this result for 
A = 0.05,0.07 and 0.1 [13]. Although, the amplitude of 
the damped oscillations decreases with increasing A, the 
rescaled spin response 5'a(a;)/5'a(0) is a universal func- 
tion of w/A* (depending on a). This is shown in the 
inset to Fig. (1) for a = 0.1. Hence, for weak dissipation, 
increasing A, for A < w^, does not destroy the damped 
oscillations. As we show below, only increasing a has 
this effect. 
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FIG. 2. 5(w) for a = 0.2 (Jy = 3.01), A = 0.001 
and different bias e/A*. A* = 0.079 A is the renormal- 
ized tunneling amplitude for the zero bias case. For the fi- 
nite bias cases, e/AJ! = ^,1,2, we find for the renormalized 
tunneling amplitude Aj/A* = 1.05,1.33,1.93, respectively, 
which is within 6% agreement of the weak coupling result [5] 
Aj, = -^Z A,.*^ -|- e^. Energies are in units of i? = ^ = 1. 

Fig. (1) for a = 0.1 and the zero bias case of Fig. (2) 
for a = 0.2 show that the damped oscillations become 
strongly suppressed with increasing a, even for A << 1. 

The oscillations disappear completely at approxi- 
mately a = 0.33 [14], well below the Toulouse point 
a = I (see below). A small bias, e << A, i.e. small 
relative to the bare tunneling amplitude, but comparable 
to the renormalized tunneling amplitude. A*, leads to a 
suppression of the inelastic peaks, and an increase in the 
tunneling frequency relative to the zero bias value. Fig. 
(2). This increase in the effective tunneling frequency. 



Ai(e), is well described by the weak coupling theory ex- 
pression [5], Ai,(e) = ^/^^TA*^. The effect of a bias, 
however, is not simply to renormalize the tunneling fre- 
quency, since we could not scale S[lli) for different e onto 
a single curve, as is the case for different tunneling am- 
plitudes and zero bias. Thus a finite bias for weak dis- 
sipation does not change the qualitative aspects of the 
dynamics. Beyond a 0.2 the inelastic peaks are negli- 
gible, even for very small values of A, and the incoherent 
part dominates (see also Fig. (4)). 

For the exactly solvable Toulouse limit, a = |, the 
numerical results for S[lli) fit very well onto the resonant 
level model result for all values of A as shown in Fig. (3). 
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FIG. 3. S[cj) for the exactly solvable Toulouse case, a = ^, 
(J|| = 1.312). A = 0.1 (circles), A = 2 4 0.1 (diamonds), 
A = 2 = 0.1 (squares). The solid lines are the exact analytic 
results for the resonant level model. 



The response consists of a single peak at uj 
^RL = t^pVrl where Vrlw^^ = 



-1 
■'km ' 



of width 
DS^(5) and 



^RL,^KM are the respective high-energy cut-offs of the 
resonant level and Kondo models [7]. Since ujkm = 2D = 
iVc, we can relate the various models by determining ujul 
from the exact resonant level model result S[lli = 0) = 



We found that 



1.125 in all cases with a 



variation of less than 0.1% [13]. 

The regime | < a < 1 corresponds to the antiferro- 
magnetic Kondo regime, for which there have been few 
reliable results for the dynamic susceptibility. The renor- 
malization group results presented here provide an essen- 
tially exact solution in this region, as can be seen from 
Table I. The single quasielastic peak in S[lli) narrows 
exponentially for a — > 1, i.e. J±^\\ — > 0, (Fig. (4)). Tun- 
neling is incoherent for all values of A. Fig. (4) shows 
the crossover from damped oscillations for weak dissipa- 
tion to incoherent relaxation for strong dissipation. The 
crossover occurs at approximately a = 0.33. From the 
scaling properties of 5'a(a;) for different A this value is 
independent of A (for A << w^). The crossover point 
separates two strong-coupling regimes, corresponding to 
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(a) damped short time oscillations for small a, and (b) in- 
coherent relaxation of the spin and Kondo behaviour for 
large a. The details of the fixed points and the behaviour 
in the ferromagnetic regime Jy < 0, — Jy > Jx, (a > 1), 
where we obtain the expected localization of the particle, 
will be discussed elsewhere. The inset to Fig. (4) shows 
that the low energy scale is given correctly by the scaling 
result, A^, for weak dissipation and by the Bethe ansatz 
expression for the Kondo temperature of the AKM [10] 
for a — > 1. 
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FIG. 4. S{uj) for A = 0.1 and different values of a. The 
inelastic peak disappears at approximately a = 0.33. The low 
energy scale, A^^°, for A = 0.1 and for a = 0.1, 0.2, . . . , 0.9 
is shown in the inset. It is defined to be the position of 
the inelastic peak, AJ! , for small a and the half-width of the 
quasielastic peak for large a. It agrees with the scaling result, 
A,., for small a, and with the Bethe ansatz Kondo tempera- 
ture for the AKM [10], Tk{J±, J\\), for a ^ 1. At a w 0.3, 
close to the crossover, there is an ambiguity in defining the 
low energy scale and some discrepancy is seen. Similarly, 
Tk{J\\, J±) deviates from the true energy scale, A^^", when 
a < 0.5 ( Tk{J\\, J±) in [10] is derived for small Jy i.e. a —> 1). 

To summarize, we have shown, that at low energies, 
0<a;^A<<a;c, the equilibrium dynamics of the dis- 
sipative two-state system is very well described by the 
AKM for dissipation < a < 1. Agreement with exact 
results was obtained in the limiting cases a — > 0, |, 1. 
The validity of the generalized Shiba relation, within er- 
rors consistent with the numerical data, was verified in 
the range < a < 1 and A << w^. 
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